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In this work we present a simple method to reconstruct the complex spectral wavefunction of 
a biphoton, and hence gain complete information about the spectral and temporal properties of a 
photon pair. The technique, which relies on quantum interference, is applicable to biphoton states 
produced with a monochromatic pump when a shift of the pump frequency produces a shift in the 
relative frequencies contributing to the biphoton. We demonstrate an example of such a situation in 
type-II parametric down-conversion (SPDC) allowing arbitrary paraxial spatial pump and detection 
modes. Moreover, our test cases demonstrate the possibility to shape the spectral wavefunction. 

This is achieved by choosing the spatial mode of the pump and of the detection modes, and takes 
advantage of spatiotemporal correlations. 


With the ability to exhibit nonclassical properties such 
as entanglement, photon pairs (or biphotons) are of 
fundamental interest in quantum optics and constitute a 
useful resource. Many proposals exploiting this resource 
require knowledge of the biphoton quantum state. The 
complex spectral wavefunction, which contains both 
amplitude and phase information, provides complete 
knowledge of the spectrotemporal state of a photon 
pair, including the ability to calculate all observables 
related to this degree of freedom and predict interference 
phenomena. Unsurprisingly then, the reconstruction of 
the full complex spectral wavefunction has received a lot 
of attention. Several approaches have been pursued over 
the last years, but each proposal entails experimental 
challenges. Interferometric methods require a high level 
of stability m- Other methods rely on nonlinear 
optical effects which are inherently inefficient at the 
low intensity levels typical of quantum light sources, 
requiring very large nonlinearities or high powers [MS]. 

In the pioneering work of Hong, Ou and Mandel 
(HOM) [7], the coherence length and time delay between 
two photons was measured using quantum interference 
on a beam splitter, circumventing the need for optical 
nonlinearities. In fact, the interference phenomenon has 
since then proven useful in a variety of applications, 
including quantum teleportation jS], quantum gates 0 
uni, linear optics quantum computation dl, Bell-state 
analysers |12j , and the measurement of the group velocity 
of light [13] , as well as of dispersion [14] . Extensions of 
the HOM approach also enable the full reconstruction 
of complex spectral wavefunctions: Chen and co-workers 
m rely on the time resolution of the detectors to directly 
measure the delay distributions, and therefore their 
method is applicable to very narrow-band biphotons. In 
contrast, Douce et al. [16] propose a scheme to measure 
the biphoton Wigner function using HOM interference by 
adding shifts of the biphoton frequencies. Yet, a practical 


implementation of such shifts is not particularly simple 

or efficient [ 13 - 

In this letter, we propose and implement a variation 
of the scheme in Ref. m that relies on the ability to 
effectively shift the relative frequency of the biphoton 
state in the generation process. As an example, our 
method allows us to measure the complex spectral 
wavefunction, and consequently also the time delay 
distribution, for type-H SPDC with a monochromatic 
pump in an arbitrary paraxial spatial mode, after 
projection of the down-converted photons into a likewise 
arbitrary paraxial spatial mode. The assumption of a 
monochromatic pump beam means that the frequencies 
of signal and idler photons are perfectly anticorrelated, 
hence reducing the problem to the determination of a 
complex-valued function of one variable. Our scheme 
is an extension of a conventional HOM type set-up 
by tuning either the temperature of the nonlinear 
crystal or the pump frequency, so that a quantum 
interference coincidence pattern is recorded as a function 
of path length difference and crystal temperature or 
pump frequency. Both pump frequency tuning and 
crystal temperature control produce the same effect: 
a frequency displacement of the wavefunction. Other 
systems, such as four-wave mixing in atomic species 
can be similarly controlled by tuning the frequencies 
of the pumps m- We show that multivariable 
quantum interference patterns can in fact be used 
to reconstruct the complex spectral mode function 
$ (H), which determines the wavefunction Idt) = 
f dfl $(H)a| -I- 0) al — O) |0). Here dl^(uj) is the 

creation operator for a photon with frequency w and 
polarisation as indicated by the subscript, and ojp is the 
pump frequency. 

We also demonstrate that in our chosen experimental 
implementation, the spectral wavefunction can be 


2 


influenced through the choice of the spatial detection 
modes, owing to spatiotemporal correlations. This leads 
to nontrivial complex spectra with marked differences 
to the standard sine function m, making their 
characterisation worthwhile. The ability to shape 
the spectral wavefunction is important for quantum 
information and communication applications, and has 
already been pursued for single photons and photon pairs 
using other approaches PUS]. 



FIG. 1. Experimental set-up: An 8 mW, monochromatic 
404.25 nm pump beam is focused to a waist of Wp = 4.3 
fim into the temperature-controlled nonlinear crystal (15 mm, 
ppKTP). The down-converted light is collimated by a lens 
after the crystal. The pump beam is discarded by a longpass 
filter. The photon pairs (s and i denoting signal and idler, 
respectively) are separated by a polarising beam splitter 
(PBS). A set of waveplates (shown as the grey line in the 
signal path) is used to maximise interference. The path length 
between the two arms differs by a controllable amount AS, 
before coupling into single-mode fibres. Alternatively, we can 
select different higher order modes with diffractive elements, 
prior to the fibre coupling. Finally, the photons pass through 
a fibre beam splitter (BS), and coincidences (&) are detected 
across two avalanche photodiodes (APDs), one for each fibre 
beam splitter output arm. 


joint function of path length difference and crystal 
temperature. Equivalently, it is possible to vary the 
pump laser frequency instead of the crystal temperature. 
We consider both cases in our theoretical analysis. 

The normalised coincidence count rates can be 
modelled as 

i?coi„c(A5,r,Wp) = t^+r^-2rh‘^Re[f {AS,T,u;p )], (!) 

where AS = {Sg — Si) is the difference between signal 
and idler path lengths, T is the crystal temperature, ojp 
the pump frequency, and t and r are the moduli of the 
transmission and reflection amplitudes of the HOM beam 
splitter, respectively. In addition, the interference term 
reads 

f{AS,T,ujp) = j dn<i>{n;T,ujp)<i>* (-0;r,Wp) 

X exp (iAS2fl/c) , (2) 

where c is the speed of light in vacuum. The conventional 
HOM dip is a slice of such a surface i?coinc(AS', T, Wp) 
along the direction, keeping the crystal temperature 
and pump frequency fixed. The measured coincidence 
counts thus involve our complex wavefunction of interest, 

$(0;T,Wp) = J y dq^dq, $fuii(qs,qnf^;r,a;p) 

xG:(q,)Gr(q.), (3) 

where $fuu(qs, qi, T’, Wp) is the wavefunction before 
projection into the spatial modes Gs(qs) and Gi(qi), 
with q being the transverse momenta ^D] . 

However, due to the nature of quantum interference, the 
wavefunction appears in the form of 


Fin, T, ojp) = d> (H; T, tCp) (-H; T, cOp) , (4) 


The experimental set-up used in the implementation 
of our reconstruction scheme is shown in Fig. 
Photon pairs are generated by pumping a periodically 
poled Potassium Titanyl Phosphate (ppKTP) crystal 
with a focused Gaussian beam in a collinear, type-H 
down-conversion configuration. They are then separated 
with a polarising beam splitter, fibre coupled, and are 
recombined on the two ports of a fibre-based beam 
splitter. As in a typical HOM experiment, the path 
length difference between the two arms can be swept. 
In addition, the temperature of the crystal is actively 
stabilised with a precision of ±20 mK and may be 
tuned to any desired value. The detection mode is 
experimentally set to optimise the singles count rates for 
the chosen pump beam waist. To create different test 
cases, we modify the usual Gaussian detection mode by 
using diffractive elements to project the down-converted 
photons into a Laguerre-Gaussian mode, or by displacing 
the nonlinear crystal along the beam’s propagation axis. 

We measure the coincidence count rate as a 


which we refer to as the symmetrised wavefunction. Since 
Fin, T, Up) is Hermitian w.r.t. n, f iAS,T,Wp) is real 
and from the coincidence rates (recall Eq. (0): 

fiAS,T,Up) = {t^ + r^-R,oinciAS,T,Up))^ . 

(5) 

We then obtain F(f2,T, Wp) by taking the Fourier 
transform of / iAS,T,ojp) with respect to AS: 


Fin,T,Up) = 


dASf iAS, T, Up) exp(-i(2H)A5/c). 


( 6 ) 

Because the symmetrisation is not isomorphic, Eq. 
(0 can in general not be inverted to retrieve the 
wavefunction from the usual HOM dip, so additional 
information is required. One possibility is to extend 
the measurements by shifting the relative frequencies of 
signal and idler. We now show that the reconstruction 
is also possible by performing a temperature or a pump 
frequency sweep. To see this, we perform a multivariate 
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Taylor expansion to leading orders of the wavevector 
z-components for the pump, signal, and idler (indicated 
by subscript to), about the values at which perfect phase 
matching takes place: at frequencies ujm = ^om, crystal 
temperature T = Tq, and transverse wavevector = 
0. From the Taylor series approximation, symmetrised 
wavefunctions at different temperatures can be related by 
shifting the frequencies, while keeping the temperature 
fixed pn| : 


‘1’ (^i Tq + AT, ojQp + Awp) To + AT, wop + Awp) 

~ $ (n + ATct + AuipCuip', TqjUIqp) 

x$* (—n + ATct + AwpC(^p;To,a;op), (7) 


where we have defined 
ct = - 
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Here, are the wavenumbers in the crystal (pump, 
signal, and idler indicated by subscripts), ojm the 
frequencies, and A the poling period of the crystal. 
The are inverse group velocities. All derivatives 
are evaluated at the reference temperature Tq and 
frequencies = Womi w G {p^s^i}. We identify Ct and 
Cujp as proportionality constants between a shift in T or 
Wp, and n. They represent the measurable sensitivity of 
the biphoton spectrum to the crystal temperature and 
pump wavelength. The complex mode function can be 
obtained as a slice through T(n, AT, Aujp) : 


$ (2ctAT + 2cuipAujp-,To,u}op) 

_ F*(—CtAT — c^pAu;p,To + AT,wgp + Awp) 

^y\F{0,TQ,uJop)\ 

( 12 ) 

Fig. i illustrates the data analysis process, when we 
choose to do a sweep in the temperature. 

In Fig. we demonstrate the use of our reconstruction 
method, with theory and experimental results for three 
test cases. The test cases all use a Gaussian pump 
beam with a beam waist of 4.3 /xm, but differ in the 
detection modes. These are (a) Gaussians with the 
crystal centered, (b) Gaussians with the crystal displaced 
by 3 mm along the propagation direction, and (c) the 
Laguerre Gaussian modes (azimuthal index, radial index) 
= (1,0), (-1,0) with the crystal centered. 

We first show the results of our theoretical analysis, 
illustrated in the left panel of Fig. 3. From our 
model of the nonlinear process and detection, we obtain 
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FIG. 2. Steps to determine 'l’(fl) from Rcoinc{AS, AT). 
Starting with Tcoinc, the real coincidence counts, we perform 
elementary operations to obtain another real function / (Eq. 
©)■ Taking a Fourier transform of / w.r.t. AS, we get to the 
complex function F (by Eq. |6f). The desired wavefunction 
is obtained by taking an appropriate slice of F (by Eq. (121). 


the expected spectral wavefunction directly (black solid 
line) [5D]. We then calculate the expected quantum 
interference pattern, i?coinc, by Eqs. Q and ([^, based on 
which the spectral wavefunction is reconstructed using 
Eqs. 0, @ and ( |I2| (red dotted line). Next we 
calculate the spectral and time delay distributions both 
from the original and the reconstructed wavefunction, 
shown in black solid and red dotted lines, respectively. 
For the experimental results, we measure the quantum 
interference pattern and reconstruct the wavefunction 
using Eqs. (§, (§ and ( |12[ ), from which the time delay 
distribution is obtained (red dotted line). 

The theoretical analysis allows us to compare the 
original and the reconstructed wavefunction, showing a 
good agreement in all cases. However, the reconstruction 
is insensitive to those quadratic and higher order phases 
as a function of fl that stem from the propagation of 
the photon pairs to the end of the crystal, or through 
any additional dispersive elements. This causes an error 
in the reconstructed phase [5D]. For any particular 
implementation, the error is limited and depends on 
the optical properties and length of the crystal, as well 
as the spectral bandwidth. The spectral bandwidth is 
influenced by the detection mode, and our choice of 
a small detection beam waist corresponds to a broad 
spectrum, which allows us to explore limitations of the 
method. The deviation in phase cannot be seen easily in 
Fig. 13 because it is comparatively small, but there is a 
difference which is quadratic in and reaches up to 0.46 
radians for the plotted section of frequencies. This, in 
turn, has a visible impact on the time delay distribution 
in (c), where a small deviation between reconstructed 
and calculated distributions is evident. We note that this 
lack of sensitivity of the reconstruction method does not 
mean that it is overall only sensitive to linear functions 
of the phase, as the phase imparted through the spatial 
projection can be arbitrary and is recovered by our 
method. 
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FIG. 3. Demonstration of the reconstruction method on three test cases. Theoretical (left) and experimental (right panel) 
results using as detection modes (a) Gaussians with the crystal centered, (b) Gaussians with the crystal displaced by 3 mm along 
the propagation direction, and (c) a pair of Laguerre Gaussian modes (1 0), (-1 0) with the crystal centered. We show, in both 
cases, (from left to right) the coincidence counts, the complex spectral wavefunctions (amplitude in arb. units), the spectral 
distributions (arb. units), and the time delay distributions (arb. units). For experiments, we omit the spectral distribution. The 
red dotted lines are the reconstruction results, while the black solid lines within the theory section are based on the simulated 
wavefunction. Our experimental conditions are simulated using To = 58 °G, so that ct = —4.8698 x 10^^ (°G-s)“^ and with a 
detection beam waist of 9.6 fim. 


The experimentally measured coincidences allow us 
to determine the complex spectral wavefunctions and 
time delay distributions of our experimental photon 
pairs. Interestingly, they also allow us to identify small 
imperfections in the experiment. For example, a slight 
off-centering of the crystal in cases (a) and (c) results 
in an asymmetry of the coincidence count map w.r.t. 
and an increased slope of the phase and mean time delay. 
In (c), we attribute differences with the theoretically 
predicted wavefunction to the fact that the radial profile 
of the theoretical detection mode is slightly different from 
the one in the experimental implementation. 

A comparison of the three rows in Fig. shows 
significant differences between the test cases. The 
Gaussian detection with the crystal centred (a) yields 
a quantum interference pattern that is symmetric w.r.t. 
AS, about a value that depends on the time delay 
acquired when signal and idler photons propagate 
through half the length of the nonlinear birefringent 
crystal. The spectrum’s departure from a sine squared 
function is highlighted by its asymmetry. It arises from 
our use of a small detection beam waist, and also results 
in the asymmetry of the quantum interference pattern 
in the AT direction. The time delay distribution has 
a symmetric peak centred at the time delay acquired 
by propagation through half of the crystal. When 
the crystal is displaced (b), the quantum interference 


pattern becomes asymmetric, the phase changes, and the 
time delay distribution shifts [2T]. Using the Laguerre 
Gaussian detection mode with the crystal centred (c) 
changes the structure of the quantum interference 
pattern markedly, even transforming the dip into a 
peak. The phase of the wavefunction is similar to 
the Gaussian case, but the spectrum has a side lobe. 
Interestingly, the time delay distribution has a dip at 
the approximate axis of symmetry, which means that 
the probability of photons arriving with their mean time 
delay is suppressed. 

In summary, we have proposed and demonstrated a 
method to reconstruct the complex spectral wavefunction 
of a biphoton, using HOM interference for type-II SPDC. 
The essence of our method lies in the fact that a change 
in temperature or pump frequency is approximately 
equivalent to a shift of the frequency for the symmetrised 
mode function that determines the quantum interference 
coincidence counts. 

A considerable advantage of the method lies in its 
simplicity, both in the experimental implementation and 
in the data analysis. Indeed, our technique is not faced 
with challenging stabilisations typical of interferometric 
measurements [3], or the need for high pump powers 
incurred by measurements that rely on nonlinear optical 
effects [ 3 . As an extension, following the results in 
Ref. [H], f {AS,T,ujp) provides the Wigner function in 
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the case of a mixed state. For the case of a pulsed 
pump where the biphoton wavefunction depends on both 
signal and idler frequencies, / (AS, T, wop) can provide 
the Wigner function in which the sum frequency variable 
has been traced out. We note that to appropriately 
manipulate the biphotons in the generation process, so 
that their wavefunction can be shifted as shown here, 
the argument of the phase matching function must be 
linear in the relative frequency fl. This includes processes 
such as collinear and noncollinear type-II SPDC with 
or without periodic poling |20| . but not type-I collinear 
degenerate down-conversion due to the group velocities 
of the two photons being equal. Moreover, post-emission 
spectral manipulation, such as the use of spectral filters 
or propagation through dispersive elements, will lead to 
a faulty reconstruction. A further limitation is that the 
reconstruction is unsuccessful at recovering the limited 
part of the quadratic and higher order phase that arises 
from the propagation of the biphoton to the end of the 
crystal [^ . 

Lastly, the freedom to choose spatial pump and 
detection modes offers some interesting possibilities. 
We have characterised the spectrotemporal properties 
of the biphoton, after projection into a spatial mode. 
Contrary to the intuitive idea that spatial degrees of 
freedom should not play a role, our results show that 
the choice of detection modes can have a pronounced 
effect on the spectrotemporal properties, in particular 
due to spatiotemporal correlations in the biphoton 
wavefunction m [2^27j . Our method works for 
arbitrary paraxial pump and projection modes, so it 
is possible to influence the detected wavefunction by 
adjusting the modes. 
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